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1 Introduction 



With this paper we continue a series of works on manifestly covariant relativistic 
statistical mechanics. In a previous paper |lj a manifestly covariant generalized rel- 
ativistic Boltzmann equation for both cases of the system of N distinguishable and 
indistinguishable events treated simultaneously has been made in the framework of a 
manifestly covariant classical and quantum mechanics [0 . These events, considered as 
the fundamental dynamical objects of the theory, move in an 8iV-dimensional phase 
space. Their motion is parametrized by a continuous Poincare invariant parameter 
r called the "historical time". The collection of events (called "concatenation" H) 
along each world line corresponds to a particle in the usual sense; e.g., the Maxwell 
conserved current is an integral over the history of the charged event Hence the 
evolution of the state of the iV-event system describes the history in space and time 
of an iV-particle system. 

A generalized relativistic Boltzmann equation was used to prove the if -theorem 
for evolution in r. In the equilibrium limit the covariant forms of the known non- 
relativistic statistical distributions were obtained. Since these distributions are the 
distributions of the 4-momenta of the events, m? = —p 2 = —p^p^ is a random variable 
in a relativistic ensemble. 

The case of distinguishable events was treated in a series of papers H,[H]> 0- hi 
[5] the simplest case of a narrow mass shell p 2 = —m 2 = —M 2 was considered, where 
M is a given fixed parameter with the dimension of mass (an intrinsic property of 
the particles), assumed to be the same for all the particles of the system. The results 
obtained in this approximation are in agreement with the well-known results of Synge 
H from an on-mass-shell relativistic kinetic theory. 

In [6] the general case < m < oo was considered and a relativistic mass distri- 
bution was obtained. It has been shown [7] that in the Galilean limit the statistical 
mechanics of distinguishable events considered in [6] passes over to the usual nonrel- 
ativistic statistical mechanics of distinguishable particles. 

In the present paper we consider the case of indistinguishable events. We obtain 
a relativistic mass distribution for this case, from the equilibrium solution of the 
manifestly covariant generalized Boltzmann equation obtained in [1], by integration 
over angular and hyperbolic angular variables. We calculate all the characteristic 
averages, obtain expressions for the energy-momentum tensor, pressure and density 
of the events, and derive the relativistic equation of state. 

We show that in the Galilean limit the theory which we discuss here passes over 
to the usual nonrelativistic statistical mechanics of indistinguishable particles. 

Finally, the criteria of validity of the theory presented in this paper are obtained. 
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2 Relativistic mass distribution 

We begin with the equilibrium relativistic distribution function obtained in [1] (we 
use the metric g^ v = (— , +, +, +) and q = q^, p = 

I e A(q)(p+p c ) 2 +B(q) 

fo(q, P ) = C{q) e _ A{q){p+pc)2 _ B{q) ^ l = C{q) lTeA{q){p+Pc)2+B{qy A(q) > 0. (1) 
The function (1) must be normalized, according to [1], as 

r r . e A{p+ Pc ) 2 +B 

n(q) = J fo(q,p)d 4 p = C(q) J d^p ^ - eA{p+Pc)2+B , (2) 

where n(q) is the total number of events per unit space-time volume in the system 
in the neighborhood of the point q. By expanding the denominator in the integrand 
into a power series and introducing hyperbolic variables 

tt A : m > 0, < 6 < vr, < if < 2tt, -oo < (3 < oo, 

we can rewrite (2) as follows: 

oo „ 

n(q) = C(g)^(±l) n+1 / d 4 pe nA{p+Pc)2+nB 

n=l 

oo , 

= C(?)X!(±l) n+1 / m 3 smh 2 (3 sin 9dmd(5d9dip 

n=l 

g— nAm 2 — 2nAm c m cosh f3— n(Am 2 — B) 

After some calculations || we obtain the normalization relation 

n(q)=C(q) 7 —J2- J ^ e ~ n( c ~ B) ^( 2 ^ Am2 c)- ( 4 ) 

ZA n=l H 

Here ^f(a, b; z) is the confluent hypergeometrical function of z (ref. [10], p. 257, sec. 6. 6). 
After integration over angular and hyperbolic angular variables we obtain from (3) 
the expression 



27T ~ (±1 )"+' r x 
l c n =l 

from which we identify the mass distribution function 



?7T °° (-\-~\ roo 

n(q) = C{q)- £ - - / dmm 2 e- n{Am ^ e~ nAm ^ K x {2nAm c m), (5) 

i\vn c i Ti Jo 



f(m) =DJ2 — m 2 e- n{Am "- B) e- nAm2 K^nAm.m), (6) 



„=i n 



3 



where 

D ~ X = S E ^f^e-"^- B ^(2,2;nAm2), (7) 
4,4 ^ n 2 v c7 

n=l 



so that /(m) is normalized according to 

f(m)dm = 1. 



In (5) and (6) K\ is the Bessel function of the third kind (imaginary argument), 
where, in general, 

K v {z) = ^e^HV{iz). 

With the help of this distribution one can obtain the local average value of an arbitrary 
function of mass <j)(m): 

(<f)(m))g = / <fi(m)f(m)dm. (9) 
J o 

Let us obtain, for example, the local average values of mass and mass squared in 
relativistic gas which represent the first two moments of the distribution (6) (factors 
m result in a change in the first argument of the confluent hypergeometric function 
and its normalization): 

_ 3tt E^ =1 {(±l)"+ 1 e-"^- fl )^(|,2;nAmg)/n§} 



2 E^° =1 {(±l) n+1 e-"^- B ^(3, 2; nAm c 2 )/n 3 } 
^E~=i{(±l) n+1 e- n ( Am '- B ^(2,2;nAm2)/n 2 }' 1 y 



More generally. 



I E~ i{(±l) n+1 e" n(Am ^ B) ^(2 + |, 2; nAm^/n 2 ^} 
! ~ * v 2 ^ ±;± V 2 ^ 2 E~=i{(±l) n+1 e- n ( Am '- B )^(2, 2; nAm 2 )/n 2 } ' 



Now, as in [5], [6], we define absolute temperature through the relation 



(12) 



2Am c = —-, (13) 
k B T' v ; 

which implies that in thermal equilibrium Am c is independent of q. Hence 

, v Svr r — E^ = i{(±l)" +1 e~ w( ^ r ~ B)l ^(|, 2; n^)/n§} 

(m) = —\ 2m c k B T ? me m . (14) 

8 V E^ =1 {(±l)" +1 e- n( ^- B) ^(2,2;n^)/n 2 } 

In the limit T — > it follows from the asymptotic formula for z — > oo [JTOj] 
(a, 6; *) ~ ^~ a [l I f]( l)" a(Q + 1} - - - (a + n - + a - b)i2 + a - b) ' ' ' (n + Q ~ h) 

n=l n-2:n 

(15) 
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that 

{m) a ^T ^w^Ty* . del 

4 E~=i{(±l) n+1 e" n(5S ^ r V" 4 } 



One can also obtain in this limit 



,E~i{(±ir +i e~ n( ^~ B) /™ 6 } 



E^ = i{(±l)" +1 e~ n( ^ _B) /« 4 } ' 
Let us now calculate the first two moments of the distribution (1): 

Jd*pp»f (q,p) = I d 4 pp \ e T A X P + fX B = E^ =1 (±l) n+1 J 
^ fd*pfo(q,p) fd*p eM T )2 +\ E^° = i(±l)" +1 /^e nA(p+Pc)2+ni? ' 

(18) 

/ n i/\ _ fd 4 pp^f (q,p) _ En=i(±V n+1 Jd 4 pp^e nA ^ 2+nB 
[P P )q Jd*pf (q,p) E^=i(±l) n+1 /^e nA(p+Pc)2+nB ' 

Define the following functions: 

F = ]T(±ir +1 / d 4 pe nA ^ 2+nB = / d 4 pf (q,p), (20) 

n=l J 
°° f±l) n+1 



n=l U J 

F 2 = £ ^ / d 4 pe nA{p+Pc) +nB . (22) 



We have 



<">• = ^ul£->®< (23) 

= ^(iisSb-^"-^-*")- < 24 > 

Calculating Fi and F 2 in explicit form and using the relations 

d n u d 



dp CjX c dm 2 ' 

d 2 d 2 d 



dpc^dpcu c c d(m 2 ) 2 dm 2 ' 

we obtain with the help of the formula (ref.[10], p. 257, sec. 6. 6) 

b ', = ~a^(a + 1, 6 + 1; z) 



» SgLi{(±l) w+1 e- B(Am °- B) ^(3,3;nAm c 2 )/n 2 } 
^ /g ~ ^ c E^=i{(±l) n+1 e-"( Am '- B )^(2,2;nAm2)/n 2 }' 1 ' 

,M\ finMnl/ S^i{(±l) w+1 e-"^- fl ^(4 > 4;nAm 2 )/n 2 } 
\P P - PcPc E^ =1 {(±l) n+1 e""( Am ?- B )^(2,2;nAm 2 )/n 2 } 
^ E^=i{(±l) n+1 e"" (Am ^ B) ^(3, 3; nAm 2 c )/n 3 } 
+ ^"E"=i{(±l) n+1 e-"( Am '- B )^(2, 2; nAm 2 )/n 2 } ' ( ^ 

As in [5], [6], we make a Lorentz transformation to the local average motion rest frame 
moving with the relative velocity 

Pc 

in order to obtain the local energy density. 
The rest frame energy is 

{E)q ~ Vl-u 2 ' 

so that 

E^°=i{(±l) Tt+1 e~" (Am '~ B) ^(3, 3; nAm c 2 )/n 2 } 
zm c E ^^ (±1)n+le _ n(Am 2_ B) ^ (2 >2;ri Am 2 )/ra 2 }" { ] 

To obtain the pressure and the density of events in our ensemble, as in [5], [6], we 
study the particle energy-momentum tensor defined by the R 4 density 

T^{ q ) = Y.jdr^5\ q - qi {r)). (28) 

i 

Using the result of [5] 

(T^(q)) q = T AV Jd 4 pf ( q ,p)^f (29) 
and the expression (26) for (p^p u ) q , we obtain 



M 

+6p^p; 



Eg 3 =i{(±l) n+1 e' Tt(Am '~ B) ^(3, 3; nAm 2 )/n 3 } 
~A En=i{(±^) n+1 e- n(Am '- B) ^(2, 2; nAm 2 )/n 2 } 



.... i/ EgLi{(±l)" +1 e~ w(Am '~ B) ^(4,4;nAmg)/n 2 } 



(30) 



E^=i{(±l) n+1 e- ra ( Am '- B ^(2, 2; nAm 2 ) /n 2 }_ 

In this expression T^y is the average passage interval in r for the events which pass 
through the small four- volume AV over the point q of R 4 . 

The formula for the stress-energy tensor of a perfect fluid has the form [5,(3.39)] 

™-=^-fr + '>8& (31) 
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where p is the pressure and p is the density of energy at q. 
According to (25), 

(32) 



hence 



and 



T Ay n(g) E^i{(±l) w+1 e-"^- B ^(3, 3; «^)/n 3 } 
P AM E^=i{(±l) n+1 e-"( Am "- B ^(2,2;nAm2)/n 2 }' 1 J 

_ QT Av n(q)m 2 c E^ =1 {(±l) ra+1 e- w ( Am '- B ^(4, 4; nAmg)/n 2 } 
P + P ~ M E^=i{(±l) n+1 e- n ( Am '- s ^(2, 2; nAm^/n 2 } ' ^ ^ 

To interpret these results, as in [5], [6], we should calculate the average (conserved) 
particle four-current, which has the microscopic form 

^(9)=E/^(9-9*W)^. (35) 

i 

Using the result of [5] 

(■/"(?)>, = Tav / d 4 p^f (q,p) (36) 
and expression (25) for we obtain 



2T Ay n(g) EgL 1 {(±l)" +1 e-^- i? ^(3, 3; nAm 2 )/n 2 } 
W W/g M Pc E^=i{(±l) ri+1 e-^ m '- i? )^(2,2;nAm2)/n 2 }" 



In the local rest frame p^ = (m c , 0), 

o 2T Ay n(g)m c E^ =1 {(±l) n+1 e-"( Am -- B ^(3, 3; nAmg)/n 2 } 

W W;<? M E~=i{(±l) n+1 e-^ Am '- B )^(2,2;nAm 2 )/n 2 }' 1 J 

Defining the density of paricles per unit space volume as 

N (q) = (J°(q)) q , (39) 
we obtain the relativistic equation of state 

No E^°=i{(±l) n+1 e- n(Am '- B) ^(3, 3; nAm 2 c )/n 3 } 



P 



2Am c E~ = i{(±l) n+1 e-"( Am "- i3 )^(3, 3; nAm 2 c )/n 2 } 



E^-i{(±l) n+1 e _n{ ^ } ^(3, 3; n^)/n 3 } 

= NoksT ni-c m • ( 4 °) 

E^ =1 {(±l)" +1 e ( ^" B) ^(3, 3; n^)/n 2 } 

In limiting case T — * the general expression for the distribution function (6) can 
be simplified. Indeed, it is seen from (17) and (27) that in this case (m 2 ) ~ (/c B T) 2 <C 



7 



(E')m c ~ m c (fcsT). Therefore, one can neglect m 2 = —p 2 in comparison with 2App c 
in the exponent of the initial distribution (1); hence, we begin with 

/o(P.i)= g M ^-B Tr ( 41 ) 

and after integration and normalization obtain 

oo /_i_i\n+l 

/°(m) =D°Y, — m 2 e- n(Am '- B) ir 1 (2nAm c m), (42) 

n=l n 

O oo Nn+1 

Cl-^E^*-" (43) 
The distribution function (42) gives 

(mr = r( | + i)r(f + 
^ r( | + i)r ( | + 2)( 2fcBr) ^^« ±1 >" +lg ""l^:'l(" w > , ,44) 

2 V 2 E^=i{(±l) n+1 e ( ^ V™ 4 } 

which coincides with the low-temperature limit of (12). Similarly, one can obtain in 
this case 

(pT = (45) 
™. E~,{(±l)" +1 e-" ( ^- B »K}' 
which coincides with the low-temperature limit of (25). 

To interpret these results, consider the initial equilibrium distribution function 
fo(q,p) in the local rest frame p^ = (m c , 0). Using (13) one obtains 

fo(q,p) = C{q) —TT—^ - - = C(q)- 



e -A(p+p c ) 2 -B ^ I ^ ^ e Am?+2Am c E+(Am%-B) ^ ]_ 

= C(q) : 5 . (46) 

^M^r[E+£- c + (^-k B T.B)}} T l 



Therefore, one can put 

M 

m c = — . 



(47) 



+ (48) 

and hence to identify (46) with the grand canonical distribution function 

1 
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obtained previously in |TTJ for the static Gibbs ensembles. In the latter formula 
fi and fix are relativistic chemical and mass potentials, respectively [11]. Hence 
Am 2 c — B = —fi/k B T, and exp{— n(Am 2 c — B)} should be identified with exp(^^) in 
all formulas containing this exponent, throughout the paper. 
For example, the relativistic equation of state (40) now reads 

Er-!{(±l) n+1 e^*(3, 3; ^fr)/n 3 } 
p = N k B T ^ mkB T>i J 4g) 

£;r =1 {(±l)" + W*( 3 , 3; ^fer)/^ 2 } 

Since at thermal equilibrium there is no dependence on q, according to [1], instead of 
fo(q,p) one can use equilibrium relativistic distribution function^ 



/o(p) = dVofap) = - ^ o - — , (50) 



normalized as 

d 4 p 



(27T)' 



Up) = n, (51) 



where n = is the event number density. 



3 Galilean limit of relativistic mass distribution 

The Galilean limit of manifestly covariant relativistic theory which we discuss here 
has been considered in a series of papers [7],[ll],|jI2|, by taking c — » oo (compared to 
all other velocities). It has been shown [7] that in the Galilean limit the relativistic 
relation between the energy E and the mass m 

E 2 = m 2 + p 2 (52) 

transforms to 

E = m + ^, (53) 

where the Galilean mass M coincides with the particle's intrinsic parameter. At the 
same time the quantity 

i] = c 2 (m — M) (54) 

may take any value, however, finite, as c —>■ oo; or, equivalently, 
m = M(l + 0(l/c 2 )). 

In case of an equilibrium relativistic ensemble of distinguishable events such a 
transformation of the relativistic relation between E and m gives rise to the usual 
nonrelativistic Maxwell-Boltzmann distribution of p 2 /2M. 

1 Clearly, all the results of the paper, since obtained for thermal equilibrium, remain valid. 
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Now we shall consider the case of indistinguishable events. We start with the 
low-temperature form (41) of initial equilibrium relativistic distribution (1) which is 
normalized as follows, 

n(q) = C(q) J d*P eAm2c _ 2A l c _ BTl - (55) 
This integral written in the local rest frame takes the form 

[ dEd 3 pjr(±l) n+1 e- 2nAmcE e- n{Am *- B) = [ dEd 3 pf^(±l) n+1 e~^ r e^ T , (56) 

^ n=l n=l 

or, in view of (53), (54) and the relation — A < 77 < A [7], 

r M+A _ JHn _ f 2_ r pi_ 1 

22 / dme I d s pe " b t^m w_ (57) 



71=1 



Af-A 



Integration on m gives 

m+a _ Jim _ k B T _hm_ nA . 

dme k B T = 2 e k B T sinh . (58) 

m-a n KbT 

Since A may take any infinitesimal valueQ, but not zero (it enters the normalization 
factor), we obtain in the limit A — > 

k B T . , nA 

sinh- £ A, 

n KbT 

and therefore 

/oo 2 r \ 
d 3 p V(±l)™ +1 e _ ^ (537 ^ +M) = 2AC(g) / (i 3 p — - 2 , 
n=l ^ eiju-ft/bBT nr 1 

(59) 

which is the usual (normalized) nonrelativistic Bose-Einstein/Fermi-Dirac distribu- 

2 

tion for e ee 

S-' = r(§)(* 1) T)V.g(^ e *. (6i) 

In the latter formulas 

p, = fi-M (62) 
is the usual nonrelativistic chemical potential. 



2 It corresponds to infinitely sharp mass shell A = c 2 | m — M | [7] . 
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A precise equation for \i = //(n, Hk,T) can be obtained from (50), (51) by per- 
forming integration: 

The corresponding equation for the low-temperature case is obtained from (41), (51), 
or directly from (63), using the asymptotic formula (15), 

(k B T) A ~ (±l) s+1 ^l. , . 

By means of the formula 

00 f±l? +1 00 (±l) s+1 00 (±z) s M 

E^ S5T = E^ = ±E i 7 L = ±«p(H * = e ^> (65) 

s=l s=l s=l 

where[] 

Li v (z) = ^2 j^i | ^ | < 1 or | z \— 1, Re is > I, (66) 
fe=i ^ 



is the so-called poly logarithm |15| , the equation (64) can be rewritten as follows, 



n = ±%^Lz 4 (±e^). (67) 

For comparison, the equation of nonrelativistic statistical mechanics for jx = jl(No, T) 
is obtained from the normalization relation of the distribution function (59), 

c/ 3 p 1 , , 

N , (68) 



(2n) 3 e (* M -a)/k B T Tl 

N being defined by (39), and is found to be 

fMk B T\ 3/2 , _g_, 
iV = ± I 2 * J Li|(±e fc sT). (69) 

4 Validity criteria 

Now we wish to consider the limits of applicability of the theory which we discussed 
in the present paper and obtain explicit criteria for its validity. 



3 For | z |> 1 the function Li v {z) is defined as the analytic continuation of this series. 
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It is clear from the general considerations that the theory is valid at medium 
and low temperatures^ still far from absolute zero, at which the degeneracy effects 
manifest themselves. The explicit criterions of its validity can be obtained as follows: 

1) Since in any system of macroscoping size the spacing between successive mo- 
mentum levels should be exceedingly small compared to the characteristic thermal 
momentum k B T, one hasQ 

P " ~ = « l (70) 

k B T LV>k B T ' 1 ' 

where is a character size of a system in /i-direction. Taking into account the four 
similar relations for \x = 0, 1, 2, 3, (70) can be rewritten as follows (in the system of 
units we use), 

N n 

T*jtiT? <<1 ' ° r (W <<1, (71) 

where L is a character size of a system, L A ~ the total four- volume, occupied 
by the system in four-dimensional space-time, and n = y^j is the event number 
density. Hence, the criterion of validity is determined by a dimensionless number S, 
conventionally called the degeneracy parameter: 

5 = « 1. (72) 

This criterion can be violated at low temperatures near absolute zero or at high event 
densities. The two possible situations constitute the so-called degenerate case and 
will be treated in the following research [13]. 

2) In order that the expansion 

1 oo 
_ \^/_|_i\n+l nA(p+p c ) 2 +nB 

p—A(p+p c ) 2 —B ^ i z^y^ 1 ) e > 

c ' n=l 

which is basic for the theory, be valid, we must have 

e A(p+p c ) 2 +B < i 

or, equivalently, 

-A(p + Pc ) 2 - B > 0. 
This relation can be rewritten in the local rest frame, using cosh/? > 1, 

Am 2 + 2Amm c cosh (3 + Am\ - B > A{m + m c ) 2 - B > 0. 

4 At high temperatures it goes over to relativistic statistical mechanics of distinguishable events 
considered previously in [6]. 

D This condition implies that in the thermodynamic limit, summation over the discrete variables 
rii can be replaced by integration over the continuous variable p, where rij is the occupation number 
in the event energy-momentum space. 
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Since the latter must be valid for any m, we deduce that 



Using the relation (48), B = + we finally obtain 

/i < 0. (73) 

This criterion is analogous to the nonrelativistic one and holds at all temperatures 
far from absolute zero, similar to the nonrelativistic case. 

3) Considering m c as a parameter, in order to determine absolute temperature for 
thermal equilibrium, one should have, in accordance with (13), 2Am c > 0. It leads to 
m c > 0, since A > 0, and, through (47), to 

Hk > 0. (74) 

The three conditions (72)- (74) constitute the set of criteria for validity of the theory 
discussed in the present paper. 

The dependencies of /i and /Ik on T are determined by equation (63) which in 
the low-temperature limit transforms to (67). The corresponding equation for the 
high-temperature case can be obtained from the normalization relation of relativistic 
statistical mechanics of distinguishable events discussed in [6], 

(2tt) 3 fi K 2(i K k B T 
by taking the limit T — > oo [6,(30)]: 

n= — {k B Tfe^. (75) 

Summarizing, we can write down the following relations, taking into account the three 
cases considered above, 

±2^3-(/c B T) 4 Li 4 (±e r BT T ), low temperatures 

™=< 7^C ( ^ T)2E - l{(±1)S+le ^^ (2 ' 2; ^ro )/s2 ^ intermediate case 

4-7T 3 jlK 



T-3—(kBT) 3 e k B T . high temperatures 



5 Concluding remarks 

We have considered an equilibrium relativistic ensemble of indistinguishable events, 
described by an equilibrium relativistic distribution function representing covariant 
generalization of the known nonrelativistic Bose-Einstein/Fermi-Dirac distributions, 
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with variable mass. We have shown that this distribution can be identified with a 
grand canonical distribution function first introduced by Horwitz, Schieve and Piron 
in their work on static Gibbs ensembles. 

We have found that an equilibrium state of such a system is characterized by 
a well-defined mass distribution, following directly from the equilibrium relativis- 
tic distribution, by integration over angular and hyprebolic angular variables. The 
normalization relation of this distribution represents an equation linking the basic 
thermodynamic properties of the system /x, n, T. 

The dependence of /i on T is explicitly defined only for low temperatures. Al- 
though the dependence of fix on T is not determined, several remarks on this point 
can be made. Miller and Suhonen studying possible applications of grand canon- 
ical distribution function of ref. [11] in high-energy physics, showed that at low tem- 
peratures one should have u K << 1, whenever at high temperatures llk >> 1- The 

2 

former relation may also justify the fact that it is possible to omit the term ^fJ? in 
comparison with (E — fi) in the exponent of the equilibrium relativistic distribution 
function in the low-temperature limit. 

Finally, we have shown that in the Galilean limit a statistical mechanics of indistin- 
guishable events discussed in the present paper passes over to the usual nonrelativistic 
statistical mechanics of indistinguishable particles. 
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